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On the basis of the general nonhnear theory of a hyperelastic material with initial stress, 
initially without consideration of the origin of the initial stress, we determine explicit 
expressions for the stress-dependent tensor of incremental elastic moduli. In considering 
three special cases of initial stress within the general framework, namely hydrostatic stress, 
uniaxial stress and planar shear stress, we then elucidate in general form the dependence 
of various elastic moduli on the initial stress. In each case the effect of initial stress on 
the wave speed of homogeneous plane waves is studied and it is shown how various special 
theories from the earlier literature fit within the general framework. We then consider 
the situation in which the initial stress is a pre-stress associated with a finite deformation 
and, in particular, we discuss the specialization to the second-order theory of elasticity and 
highlight connections between several classical approaches to the topic, again with special 
reference to the influence of higher-order terms on the speed of homogeneous plane waves. 
Some discrepancies arising in the earlier literature are noted. 



Keywords: elastic moduli, isotropic stress, initial stress, invariants, plane waves 



1 



1 Introduction 



Residual stresses in solids, i.e. stresses that are present in the absence of load (body forces 
and surface tractions) can have a very significant effect on the mechanical behaviour of the 
structures in which they reside. This is the case for materials as diverse as hard engineering 
and geological materials and soft solids such as gels and biological tissues. Equally, stresses 
that are generated due to applied loads, associated with finite deformations and commonly 
referred to as pre-stresses, have a significant effect on subsequent material response, leading 
to very different results compared with the situation in which there is no applied load. 
Residual stresses and pre-stresses are examples of initial stresses but are different in nature 
in the sense that residual stresses are necessarily non-homogeneous while pre-stresses may 
be either homogeneous or non-homogeneous. Also, pre-stresses are usually associated 
with an elastic pre-deformation, while residual stresses can result from processes such as 
manufacturing, plastic deformation, growth and remodelling, for example. In either case 
it is important to be able to analyze the effect of the initial stress on the properties of the 
material and on its mechanical response. In this paper we are concerned primarily with 
the effect of initial stress on the material properties in general, and on elastic 'constants' in 
some specific cases, as well as its effect on the speeds of propagation of homogeneous plane 
waves. The study is conducted with a view to the non-destructive evaluation of initially 
stressed solids, which are ubiquitous in Nature and Engineering. For this purpose the 
initial stress is included in the constitutive description of the material without reference 
(initially) to any finite deformation with which it may be associated, and we emphasize 
that in general the initial stress appears in a highly nonlinear form. 

The origi ns of an e lastic ity theory incl uding initial stres s can be traced back as far as to 



the works of ICauchvl (118291). accord i ng to iTruesdelll f 



966r). Notable ea rl y contributors t o 



the subje c t incl ude iPoincarg (118921 ) , iHadamardl ( 119031 ) , iRavleighl ( 119061 ) , iBrillouinI (119251 ) , 
and iLovd (Il927h . 



In the context of the modern linear theory of elasticit y, the effect of initial stress was 
first examined in the work of Biq t for static problems (iBiotl. 119391) an d also for wave 
propagation problems ( IBiotl Il940al ): see also his monograph (iBiotl . Il965f ). Here, for later 
reference, we write the components of Biot's elasticity tensor as Bpigj with respect to a 
Cartesian coordinate system (xi, 0:2, X3). In general these enjoy the minor symmetries 



B = B = B 

i~^piqj '^ipqj '^pijqy 



but when there exists a strain-energy function there is also the major symmetry connection 



B - B 

'-'piqj '-^qjpi 



"ipijq 



'^jq ' ipy 



:i.2) 



where Tj,- are the components of the initial Cauchy stress and 6ij is the Kronecker delta. 



For infinitesimal strains Cij, with Cij 



[u 



+ Uj^i)/2, where 



u 



dui/dxj and Ui are the 



components of the displacement vector, the associated Cauchy stress, additional to the 

Biot left the dependence of Bpigj on the initial 



initial stress, is given by cr^ 
stress unspecified for most of his theoretical development and he was not concerned with 



pi 



BpiqjUj^q. 
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the source of the initial stress. In this sense his theory may be referred to as the general 
linear theory of elasticity with initial stress. Biot was more specific in particular cases, 
where he considered both isotropic and planar orthotropic specializations and an initial 
stress due to hydrostatic pressure, uniaxial compression, or gravity. In particular, for his 
isotropic model Bpiqj may be written in the form 



piqj 



:i.3) 



where /io and Aq are the notations that we shall use in this paper for the classical Lame 
moduli of linear isotropic elasticity. This form did not appear explicitly in Biot's work, 
as far as the authors are aware, but may be deduced from the plane strain expressions in 
equation (8.31e) of Biot's book ( iBiotl . Il965l ). Biot did in fact acknowledge that in general 
the elastic response in the presence of initial stress is not isotropic, but he adopted an 
isotropic constitutive description for simplicity. Note that (11. 3p satisfies the conditions 

The works of Biot have formed the basis for many contributions to the literature, 
particularly relating to wave propagation problems in t he geophysic al context, which was 
the original context in which the theory was developed by lBiotl fll940al ). Other contributions 
to the analysis of initial stress, more specifica lly residual s tr ess, haye app e ared in a serie s 



of pap e rs by Hoger and co-workers, incl u ding iHogerl (|1985[) , iHogerl (Il986l ) , iHogerl (jl993al ) , 
Hogerl (jl993bl ) , iJohnson fc Hogerl (119931 ) , iHogerl (119961 ) , some of which are con cerned with 
the combined effec t of f inite deformation and residual stress, and in papers by iMan fc Lu 



(119871 ) ■ iManl (119981 ) and ISaravananl (120081 ). for example. More recently, a gen eral theory of 



nonlin ear h yperelasticity fo r an in itially stressed solid has been developed bv lShams et al. 



(120111 ) and lOgden fc Singh! (120111). the latter being focused on fibre reinforced materials. 
We also mention the paper by iBazantI (jl97ll ). who detailed connections between several 
earlier formulations of linear elasticity with initial stress and their implications for the 
stability of elastic bodies. 

We distinguish between the above approach and that concerned with the effect of pre- 
stress that is associated with a finite deformation, a subject that has attracted many 
contributions, primarily concerned with the effect of the finite deformation on the propa- 
gation of small amplitude elastic waves and associated static bifurcation problems. This is 
commonly know as the theory of incremental (or small) deformations suyerimyos ed on a 
finite defor mation . We s hall not discuss this extensive topic in detail but refer to lOgden 
(119841 ) and lOgdenI (120071 ). for example, for pointers to the literature. 

A special case of the theory of finite elastic deformations in which the strains are small 
but the linear theory is no longer adequate is sometimes referred to as second-order elas- 
ticity. In this theory the strain-energy function is expanded to the third order in some 
suitable measure of strain and the stress is second order in the strain. Most commonly it is 
the Green (or Green-Lagrange) strain tensor that is used and for an isotropic material the 
strain energy is expres sed in ter r as of invariants of the strain. The first such contribution 
appears to be that of Brillouir] ( 1925 ). although this has not always been acknowledge d 
appropriately. Equivalent formulations were developed later by lLandau &: Rumerl (119371 ). 
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MurnaghanI (119371 ). iBiotI (Il940bh. iToupin fc BernsteinI fll96lh andlHaves fc RivlinI (Il961f l: 
see also the books by Brillouin ( 1946| ). BiotI ( 1965 ) and Landau fc Lifshitz J 1986 ). There is 
also a detailed historical discussion in the volume bv lTruesdell fc Nolll fll965l ). in particular 
in section 66 therein. An important objective within these works was to determine the 
correction to the speeds of waves due to the nonlinear terms in the stress. In particular, 
various formulas were found that highlighted the effect of an initial hydrostatic pressure (as- 
sociated with a pure dilatation) or uniaxial compression on the speeds of longitudinal and 
transverse waves. Of other works dealing wi t h wav e speeds based on this weakly nonlinear 
theory we mention those of iHughes &: Kellvi (119531 ). who obtained experimental results for 
longitudinal and transverse wave speeds in polystyrene, iron, and pyrex for separate initial 
stresses corresponding to hydrostatic pressure and simple compress ion and related their 
result s to formulas based on Murnaghan's second-order theory, and iThurston fc Brugger 
(Il964j ). who obtain ed expression s for the s e cond- order corrections to wave speeds in cubic 
crystals. Papers by iBirchI (119381 ) and iTangi (119671 ) also made use of the second-order theory 
but failed to include the second-order constants in their expressions for the wave speeds 
(note that the second-order constants are sometimes referred to as third-order constants). 
We discuss this and other deficiencies of the latter paper in Sections 6.1, 6.3 and 7.3. 

For several researchers it seemed important to show that all elastic materials would 
behave similarly under an initial stress. Hence Biot's incremental moduli ( II. 3p change 
from one material to another with changes in the values of Aq and fio, but t he effect 
of the pre-stress r remains the same across all solids. Similarly, according to [Lazarus 



(ll949h . lLovel (Il927h showed that under external pressure P, the (linear) elastic constants 
C44 and C12 (Voigt notation) of a certain class of solid are linked approximately by the 
'Cauchy relation' C44 = cu — 2P. This type of behaviour under pre-stress would in turn 
lead to a corresponding effect on the shift in speed experienced by an acoustic wave (for 
instance a wave should always travel faster in a pressurized isotropic solid than in its 
unstressed counterpart). However intuitive this expectation might be, it is not supported 
by experimental facts, as confirmed by the data shown in Tables [T] and [21 which show that 
the wave speed can increase or decrease with pressure, depending on the material. 

The purpose of the present paper is to draw together and highlight connections between 
some of the his torical results within a common and fairly general framework based on the 
development of [Shams et all (120 111 ) concerned with the constitutive law of a hyperelastic 
material with initial stress. In particular, we examine how the elasticity tensor depends 
nonlinearly on initial stress, with emphasis on the important special cases of hydrostatic 
initial stress, uniaxial initial stress, and initial shear stress. 

In Section [2] we summarize the basic equations for an elastic material for which the 
strain-energy function depends on an initial stress as well as the deformation from the 
initially stressed reference configuration. In particular, we express the strain energy as a 
function of combined invariants of the initial stress and deformation and give expressions 
for the nominal and Cauchy stress tensors. Next, in Section [3l we derive the equations of 
motion for small displacements from a homogeneously deformed configuration when the 
initial stress is uniform, which leads to the need for expressions for the elasticity tensor. 
Such expressions are given in Section [H but now specialized to the undeformed (but initially 
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Table 1: Initial variation of the squared wave s peeds for several solids under hydrostatic 



pressure P, as collected by Johnson et al. (199^) : pr is the mass density at P = 0; vt and 



vl are the speeds of the transverse and longitudinal waves, respectively 




Alumina 
Aluminum 
Armco-Iron 
Fused silica 
Gold 

Magnesium 

Molybdenum 

Nickel-steel 

Niobium 

PMMA 

Polystryrene 

Pyrex 

Steel (Hecla) 
Tungsten 



1.12 
2.92 
5.7 

-1.42 
0.90 
1.47 
1.05 
1.55 
0.29 
3.0 
1.57 

-2.84 
1.46 
0.70 



4.46 
12.4 
9.3 

-4.32 
6.4 
6.89 
3.48 
2.84 
6.18 
15.0 
11.6 

-8.6 
7.45 
4.58 



Table 2: Initial variation of the speeds of transverse waves for several solids under uniaxial 
strain e: Vi2, f2i are the speeds waves travelling in the direction of tension, and orthogonal 
to the direction of tension, respectively; and ar e their values at e = 0. The first row 



of data is from experiments on a sample of rail steel ^Egle & Brau. \l97d): the other ro ws 



are from experiments on soft solids with different compositions ( Gennisson et al. . 200'i) 



Solid 




=0 



Rail steel 1 
Agar- Gelatine 1 
Agar- Gelatine 4 
Polyvinyl acetate 1 
Polyvinyl acetate 3 



-0.15 
0.84 
4.53 
0.71 
1.69 



-1.50 
-0.92 

2.26 
-0.47 

0.35 
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stressed) reference configuration. Several special cases are considered in which the initial 
stress is either purely isotropic, uniaxial, or a planar shear stress. These formulas are 
then used in Section [5] to define relevant elastic moduli that depend of the initial stress 
in question and to make that dependence explicit. The moduli include stress-dependent 
Lame moduli in the case of isotropic initial stress, Poisson's ratios and Young's moduli for 
uniaxial initial stress, and planar Poisson's ratios and Young's moduli for planar initial 
shear stress. 

In Section inithe results are applied to infinitesimal wave propagation and related to some 
known results as special cases. Section [7] then considers the deformation of an isotropic elas- 
tic material from a stress-free reference configuration in order to make contact with results 
in the preceding sections by considering, in particular, a pure dilatation and a deformation 
corresponding to simple tension. We then focus on the specialization to second-order elas- 
ticity, in which the strain-energy is approximated as a third-order expansion in the Green 
strain tensor in order to highlight connections with the theory of elasticity with initial 
stress herein a nd to draw togeth er various contributions from the literature that date back 
to the work of iBrillouinI (jl925[ ). with particular reference to expressions for longitudinal 
and transverse wave speeds. 



2 Elasticity in the presence of initial stress 

We consider an elastic body that is subject to an initial (Cauchy) stress r in some reference 
configuration, which we denote by Br- In the absence of intrinsic couple stresses r is 
symmetric. Let X be the position vector of a material point in Br and let Grad and Div 
denote the gradient and divergence operators with respect to X. If there are no body 
forces then r must satisfy the equilibrium equation Divr = 0. For the most part we shall 
not be concerned with how the initial stress arises, but in Section [7] we shall relate r to an 
underlying finite deformation and r is then considered to be a pre-stress, which requires 
appropriate tractions on the boundary dBr of Br- However, if the traction on the boundary 
dBr of Br vanishes point wise then r i s referred to a s a residual stress. A residual stress is 



necessarily non-uniform flHogerl . Il985l : lOgdenl . l2003l ) and in general is not associated with 



a deformation from a stress-free configuration. 

Next, we consider the body to be subject to a finite elastic deformation from Br into a 
new configuration B with boundary dB so that the material point X takes up the position 
X in i3 given by x = x(X), where the vector function x defines the deformation for 
X G Br- The deformation x is required to be a bijection and to possess appropriate 
regularity properties, which we do not need to specify here. The deformation gradient 
tensor, denoted F, is defined by F = Gradx, from which are formed the left and right 
Cauchy-Green deformation tensors, defined by 

B = FF^, C = FTF, (2.1) 

respectively. 
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We denote by a the Cauchy stress tensor in the configuration B and we suppose that 
there are no body forces, so that the equihbrium equation diver = holds. We shall also 
make use of the nominal and the second Piola-Kirchhoff stress tensors, denoted S and T, 
respectively, which are related to cr and each other by 

S = JF^V = TF^, T = JF^VF"^ = SF-^, (2.2) 

where J = detF > 0. The nominal stress S satisfies the equilibrium equation 

DivS = 0. (2.3) 

In the absence of intrinsic couples, cr, and hence T, is symmetric, while in general S is not 
symmetric and satisfies 

FS = S^F^. (2.4) 

We now consider the elastic properties of the material to be characterized in terms of 
a strain-energy function, defined per unit volume in B^^ which we denote by W . We write 

W = W{^,t) (2.5) 

to refiect the dependence not only on the deformation gradient but also on the initial stress. 
By objectivity W depends on F only through C = F"^F, but it is convenient to retain the 
functional dependence indicated in (12. 5p . In general the presence of the initial stress will 
generate anisotropy in the material response relative to and r has a role similar to that 
of a structure tensor associated with a preferred direction in an anisotropic material. An 
exception to this arises if r is an isotropic stress. If r is non-uniform then the material 
is necessarily inhomogeneous, but if r is independent of X the material is homogeneous 
unless its properties depend separately on X. In the present paper we shall consider r to 
be uniform and the material to be homogeneous. 

We shall not consider internal constraints such as incompressibility, in which case the 
nominal stress is given by 

dW 

S = ^(F,r), (2.6) 

and the Cauchy and second Piola-Kirchhoff stresses can be obtained from (12. 2p . When 
evaluated in Br^ (12.61) reduces to 

dW 

r=3|r(I.r). (2.7) 

where I is the identity tensor. Equation (12. 7p imposes a restriction on the admissible forms 
of strain-energy function for an initially-stressed elastic material. 



2.1 Invariant representation of the strain energy and stresses 

The strain-energy function W depends on C and r, both of which are independent of 
rotations in the deformed configuration B. Thus, W is automatically objective. If the 
material possesses no intrinsic anisotropy relative to B^-, so that it would be isotropic 
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relative to Br in the absence of initial stress, then W is an isotropic function of C and r, 
i.e. 

ly (QCQ^, QtQ^) = W{C, r) for all orthogonal Q, (2.8) 

and it can be expressed as a function of the invariants of C and r. We list a possible (and 
complete) set of independent invariants as 

/i=trC, J2 = |[(trC)2-tr(C2)], /g = det C, (2.9) 
trr, tr(T^'' ^-^^-3^ 
/6 = tr(rC) 



/7 = tr(TC2) 



J8 = tr(T^C), /9 = tr(T^C^ 



(2.10) 
(2.11) 



where we have used the standard notatio n /i , 1^1 for the p rincipal invariants of C and 
followed the notation /g, . . . , /g adopted bv lShams et al\ fl201l[ l for the combined invariants 
of C and r. In the reference configuration Br these reduce to 

/i = /2 = 3, /3 = 1, /6 = /7 = trr, /8 = /9 = tr(r2). (2.12) 

For full discussion of invariants of tensors we refer to Spencer ( 197l[ ) and Zheng ( 19941 ). 
Here there are 10 independent invariants of C and r in general, a number that may be 
reduced in a two-dimensional specialization or for specific simple deformations and/or 
initial stresses. We have not attributed notations to the invariants fl2.10p since they are 
independent of the deformation and do not contribute explicitly to expressions for the 
stresses. However, W may depend on f l2.10p implicitly, but we do not list them in the 
functional dependence and we write W = W{Ii, I2, h, h, Ij, h, h), retaining the notation 
W, which is used severally for the (F,t), (C,t) and (Ji, I2, I3, Iq, I7, Is, I9) arguments. 
From (12. 6 1) the nominal stress tensor may be expanded as 



dW 



(2.13) 



where we have used the shorthand notation Wi = dW/dIi, i G X, and X is the index set 
{1,2,3,6,7,8,9}. We emphasize that although their derivatives with respect to F vanish 
the invariants fl2.10p are included implicitly in the functional depende nce of W in general . 
The required expressions for dli/dF were given in Appendix A of (jShams et all 1201 ll ) 
and are not repeated here explicitly but used implicitly. The resulting expression for the 
Cauchy stress cr is obtained from (12. 2p in the form 



JcT = 2WiB + 2W2{IiB - B^) + 2/31^3! + 2Wei: 

+ 2Wy{1:B + BS) + 2W8SB-^S + 2Wg{i:B-^T;B + BSB^^S) 



(2.14) 



wherein we have introduced the notation S = FrF^ for the push forward of r from Br to 
B, and B = FF""" is the left Cauchy-Green tensor. 

When fl2.14p is evaluated in the reference configuration it reduces to 



2{Wi + 2W2 + Ws)I + 2{We + 2W7)t + 2{W8 + 2Wg)T^ 



(2.15) 
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where Wi, i G X, are evaluated for the invariants given b y (12.121). Equ a tion (12.15]) is the 
speciahzation of (12. 7p for the invariant form of W. As in IShams et all (120111 ) we deduce 
that 

1^1 + 2^2 + 1^3 = 0, 2{We + 2W7) = l, 1^8 + 21^9 = (2.16) 

in Br- 



3 Incremental motions 

We now consider the static finite deformation x = x(X) that defines the equihbrium 
configuration B to be followed by a superimposed incremental motion x(X, t), where t is 
time. A superposed dot signifies an incremental quantity and the resulting incremental 
equations are linearized in the increments, which are considered appropriately 'small'. 
Thus, X represents a small displacement from x. We shall also write the displacement 
in Eulerian form as u = u(x, t), noting that x(X, t) = u(x(X),t). The corresponding 
increments in the deformation gradient F and J = det F are given by the standard formulas 

F = LF, j=JtrL, (3.1) 

where L = gradu is the displacement gradient. 

The (linearized) incremental nominal stress is written 

S = AF, (3.2) 

where 

is the fourth-order elasticity tensor and, in component form, ^F = Aai/3jFjp defines the 
product used in ( 13. 2p . The usual summation convention for repeated indices is adopted here 
and henceforth. For full discussion of the theory of incrementa l defor mations and motions 



superimposed on a finite deformation we refer to lOgdenI (1 19841 . 120071 ) . for example. 

By taking the increments of the connections Jcr = FS and S = TF""" from ( 12. 2 p we 
obtain, after a little rearrangement. 

So = J^^FS = 6- + (trL)<T-L(7, (3.4) 
to = J^^FTF^ = So - o-L^, (3.5) 

wherein the notations So and To are defined. These are the updated forms of S and T, 
respectively, referred to deformed configuration, and otherwise know as their 'push forward' 
forms. The corresponding push forward Aq of the elasticity tensor is such that Sq = A^h. 
It then follows from the symmetry of cr and its increment that Sq + Lcr is symmetric, and 
hence 

AoL + L(T = {AoLf + (tL^. (3.6) 
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In component form ^ is related to A. via 



Note that as well as possessing the major symmetry Aopiqj = Aoqjpi, which follows from 
fl3.3p and {\3.7\i . ^ has the property 

which can be deduced from (13. 6p . 

We assume that there are no body forces. Then, the incremental motion is governed 
by the equation 

DivS = PrX,ti, (3.9) 

where pr is the mass density in Br and a subscript t following a comma signifies the material 
time derivative, i.e. the time derivative at fixed X, so that x.^ = u t is the particle velocity 
and X = u the acceleration. 

Equation (13.91) may be updated (i.e. pushed forward) to the configuration B by writing 
it in terms of So and u, which yields 

divSo = pntt, (3.10) 
where p = PrJ~^ is the mass density in B, or, in (Cartesian) component form, 

{AopiqjUj^q)^p = pUi^tt- (3-11) 

Henceforth in this paper we assume that the initial stress r, the underlying deformation 
F and the material properties are homogeneous, so that .4. and ^ are independent of X. 
The equation of motion (13.111) then becomes 

AopiqjUj^pq = pUi^tt- 

(3.12) 

This will be used in Section [6] in discussion of the propagation of homogeneous plane waves, 
but before proceeding to that analysis we obtain explicit expressions for the dependence of 
the components of the elasticity tensor and of various elastic moduli on the initial stress 
based on the invariants of the right Cauchy-Green deformation tensor C and the initial 
stress tensor r discussed in Section 12. 1[ 



4 Expressions for the elasticity tensor 

The elasticity tensor ,4. in (13. 3p may be expanded in terms of invariants as 

id i,j £l 
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where Wij = d'^W/dlidlj, i,j G X, and X is again the index set {1, 2, 3, 6, 7, 8 , 9|. 



The detailed (lengthy) expressions for the components of ^ were given by lShams et al. 



(120 111 ) for a general deformed configuration based on expressions for the second derivatives 
of the invariants, which were given in Appendix A of the latter paper. Here we require only 
their speciali z ation to the (undeformed) reference configuration [B — )■ Br), which, following 



Shams et al\ (1201 ll ). yields 

•^Opiqj (-^li^^ij^pq ~l~ ^iq^jp) ~l~ (-^2^ip^jq ~l~ ^ij'^pq ~l~ f^li^^ij'^pq ~l~ ^pq^ij ~^ ^iq^jp ~l~ ^jp^iq) 
~l~ /^2(^^ip'^jq ~l~ ^jq^ip) ~l~ f^Z'^ip'^jq ~l~ lli^^ij'^pk'^kq ~l~ ^pq'^ik'^kj ~l~ ^iq'^jk'^kp ~l~ ^ jp'^ik'^kq) 

+ l2{SipTjkTkq + SjqTikTkp) + 'Jsi^ipTjkTkq + TjgTikTkp) + 'y^Tik'^kpTjlTlq, (4.2) 

where the a's, /3's, and 7's are defined by 

ai = 2{Wi + W2), (X2 = 4:{Wn + ^Wu + 4^22 + 21^13 + 41^23 + W^ss -W^- W2), 
l3i = 2W7, /32 = 4:{Wie + 2Wi7 + 2W26 + m27 + W3e + 2Ws7), 71 = 2^9, 

= 4(1^66 + 41^67 + 41^77), 72 = 4(1^18 + 21^19 + 2W^28 + 41^29 + W^38 + 21^39), 
73 = 4(1^68 + 21^69 + 21^78 + 41^79), 74 = ^{Wss + 41^89 + 41^99), (4-3) 

all derivatives Wi and Wij being evaluated in the reference configuration and use having 
been made of the connections fl2.16p . Note that in general the expressions ( 14. 3 p may 
depend on the invariants trr, tr(T^) and tr(T^). There are nine (r-dependent) material 
parameters in the above, j ust as there are nine constants for an orthotropic linearly elastic 
material (see, for example, Tingl . 1996[ ). but additionally here the components tij of r are 
present separately. 

Note that when referred to axes that coincide with the principal axes of r, the only 
non-zero components of (14. 2 p are given by 

Aoiiii = 2ai + a2 + (1 + 4/3i + 2/32)n + ^st^ + 2(271 + 72)r?^ + 273rf + 74r/, 

•Aoiijj = a2 + P2{ri + Tj) + P^TiTj + 72(r^ + r|) + 73(ri + Tj)TiTj + 74r?^r|, 

•Aoijij = ai + Ti + Pi{Ti + Tj) + 7i(r^ + r|) = Aoijji + r^, (4.4) 

where there are no sums on repeated indices, i ^ j, and the principal 

values of r (in general, there are 15 non-zero components of ,4o in total, dependent on 
nine material parameters and three principal stresses). In the linear specialization of the 
above only the parameters ai, 0^2, /32 are retained along with ri,r2,r3, although we 
should strictly expand ai and 02 as linear functions of r. Then there are 9 constants 
involved, specifically a;i(0), 02(0), /3i(0), /32(0), a;i^j(0), a2,j(0), i = 1,2,3, and the three 
principal initial stresses ri, T2, t^, where = (0,0,0) is the value of (Ti,r2,r3) for zero 
initial stress and signifies differentiation with respect to Tj, i = 1,2,3. Note that since 
the coefficients in (14. 4p are symmetric functions of {ti,T2,T3) the constants ai,j(0) and 
^2,1(0) are independent of i. 

We now specialize the above to consider three specific forms of r, corresponding to 
isotropic initial stress, uniaxial initial stress and planar shear initial stress. 
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4.1 Isotropic initial stress 

Suppose that r is isotropic and write r = rl, where I is again the identity tensor and 
r > (< 0) corresponds to hydrostatic tension (pressure). Then equation (14 .21) reduces to 
the compact form 

Aopiqj = rSijSpq + a{T){5ij5pq + 5iq5pj) + (3{T)5pi5qj, (4.5) 

where the notations Q;(r) and /3(r) are defined by 

a(r) = ai(r) + 2r/3i(r) + 2Ai(r), (4.6) 
/3(r) = a2{T) + 2r/32(r) + r2[/33(r) + 272(r)] + 2rS(^) + ^Sl^), (4.7) 

and we note that, by virtue of the speciahzations trr = 3r, tr(T^) = 3r^,tr(T^) = 3r^, 
the coefficients ai, ...,74 are now (in general) functions of the single parameter r, which is 
indicated above by inclusion of the argument r. 



4.2 Uniaxial initial stress 

Here we take r = ra ® a, where a is a fixed unit vector (the direction of the uniaxial 
stress, which is a tension for r > and a compressive stress for r < 0). In this case the 
components of .4o may be expressed in the form 

Aopigj = ai{6ij6pq + 6iq6pj) + a2Spi6qj + r(/3i + T'^i){5ijapaq + Spqa^aj + Sigapaj + Spjaittq) 
+ rSijUpaq + r(/32 + r-f2){Spiaqaj + dq^apai) + t'^{(3^ + 2x73 + r^74)apaiagaj. (4.8) 

Again the coefficients cti, 74 depend on r in general but for the sake of brevity this has 
been left implicit here. Note that, in addition to r, there are five separate (combinations 
of) parameters, namely ai, 0^2, Pi + T71, [32 + T72, /Ss + r73 + r^74 and we recall that in 
classical transversely isotropic linear elasticity there are five material constants (see, for 
example, iTingi . Il996l ). 

Without loss of generality we may take a to coincide with the axis ei. Then the 
components of are listed as 

Ann = 2ai + 02 + (1 + 4/3i + 2/32)r + (471 + 272 + + 273x2 + ^^t\ (4.9) 

Am = 2ai + ^2, Aoiijj = a2, Anii = "2 + (32r + 72^^^ i,j e {2, 3}, i ^ j, (4.10) 

Aim = Aiiii = Aini = «i + + 7ir^ z G {2, 3}, (4.11) 

Ami = ai + (1 + /3i)r + 7ir^ Aoijij = Aoijji = ai, i, j e {2,3},i ^ j, (4.12) 

for later reference. Note that in the linear specialization there remain only seven indepen- 
dent constants, namely q;i(0), 02(0), /3i(0), /32(0), a'liO), 02(0), with argument r = 0, and 
r, where the prime indicates differentiation with respect to r. 



12 



4.3 Planar shear initial stress 

Consider planar shear stress in the (xi,X2) plane of the form r = r(ei (X> 62 + e2 (g) ei). 



Then the only non-zero components of .Aq are written as 

Ann = A2222 = 2ai + aa + 2(271 + 72)r^ + 74r'^, A3333 = 2ai + aa, (4.13) 

An22 = "2 + 272x2 + 74r^, An33 = A2233 = "2, (4.14) 

Aoijij = Aoijji = ai + {(33 + 27i)r^ i,j e {1, 2}, i 7^ j, (4.15) 

Aoi3i3 = Ao3i3i = Ao3ii3 = + 7l^^ ^ ^ {1, 2}, (4.16) 

Aoiiij = (2/3i + /32)r + 73r^ Aim = Amj + t, j G {1, 2}, i j. (4.17) 



In the linear specialization there are now five constants, ai(0), 0^2 (0), 2/3i(0) + /32(0), a;'i(0), 
0^2 (0), in addition to r. 

5 Dependence of elastic moduli on initial stress 
5.1 Isotropic initial stress 

In the connection (13.41) we now specialize the Cauchy stress cr to the initial stress r, so 
that 

tr = So + Lr - (trL)T. (5.1) 
It follows on use of (14. 5 p that, for an isotropic initial stress r = rl. 

So = AL = Tl7 + a(r)(L + L^) + /3(r)(trL)I, (5.2) 

and hence that 

CT = [/3(r) - r](trL)I + [«(r) + r](L + L^). (5.3) 

Since the initial stress is purely isotropic we can therefore identify the stress- dependent 
Lame moduli, which we denote as A(r) and /^(t), as 

A(r) = (5{t) -r = a2(r) + 2r/32(r) + A^r) + 272(r)] + 2r=^73(r) + t%{t) - r, (5.4) 
/i(r) = a(r) + r = «i(r) + 2r/3i(r) + 2r27i(r) + r. (5.5) 

For incremental simple shear we may, without loss of generality, restrict attention to 
the {xi,X2) plane. If the shear is in the xi direction with amount of shear Lij = Li26ii62j 
then the corresponding incremental nominal and Cauchy stress components are equal and 
given by 5*021 = 0-12 = A2i2i-Z^i2 = /^(t)-Z^i2- Similarly, for incremental simple shear in the 
X2 direction with Lij = L2iS2i6ij we have 6*012 = ctu = Ai2i2-f'2i = 1^{t)L2i, and from 
dlSl) and ([53]), A2121 = A1212 = At(r). 

For incremental pure dilatation with Lu = L22 = L33 = (trL)/3 we obtain 

ti& = [3A(r) + 2;u(r)]trL, (5.6) 
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and this enables us to define the stress- dependent hulk modulus h{t), analogously to the 
classical formula, as 

^{r) = X{T) + ^f,{r). (5.7) 

Note that from (15. ip we obtain trSo = tr6- + 2rtrL, but use of trSo instead of tro- does 
not give the correct form of the bulk modulus. Also, by using (13. 5 p with cr = t = tI and 
(15. 2 p we obtain 

to = «(r)(L + LT) + /3(r)(trL)I. (5.8) 

Thus, it is clear that because of the dependence on the initial stress, different choices of 
stress measure lead to different possible definitions of the stress-dependent elastic moduli. 

The correct definitions for the stress dependent Lame moduli are (15. 4p and (15. 5p . If the 
second Piola-Kirchhoff stress is used instead, then /i(r) and A(r) would be rep laced by ajr ) 
and /3(r), respectively. This was effectively what was done in the paper bv iTaneJ (Il967f l. 



although he worked in terms of Young's modulus and Poisson's ratio. This identification 
of the Lame moduli leads to erroneous results for the speeds of homogeneous plane waves, 
as we shall show in Section |6l 

If there is no initial stress (reduction to the classical case) we denote the classical moduli 
by XQ,fiQ,KQ, so that Aq = A(0) = 02(0), fiQ = fi{0) = ai(0), and kq = k(0). If the initial 
stress is small in magnitude then we may linearize the expressions (15. 4p . (15. 5p and (15. 7p to 
obtain 

A(r) ^ Ao + [a^(0) + 2/^2(0) - l]r, (5.9) 
/x(r) ~ /io + K(0) + 2/3i(0) + l]r, (5.10) 
k{t) ~ /to + [2a'i(0) + 3a^(0) + 4/3i(0) + 6/32(0) - l]r/3. (5.11) 



5.2 Uniaxial initial stress 

When the initial stress is uniaxial the subsequent incremental response is transversely 
isotropic in nature. Then it is appropriate to work in terms of Young's moduli and Pois- 
son's ratios. In order to determine these we need to examine both triaxial incremental 
deformations without shear and separate incremental shear deformations. Consider first 
the normal components of L, written Ln, L22, -^^33, and take ei to be the direction of the 
uniaxial initial stress. Then, bearing in mind the symmetry in (I4.10p . the corresponding 
incremental nominal stresses are 

>S'oil = ^oiiiiLii + ^0ll22(-^22 + -^33)) (5.12) 
>S'o22 = -4.01122-^11 + -4.02222-^22 + -4o2233-^33) (5.13) 
>S'o33 = -4oil22-^^ll + -4o2233-^22 + -4o2222-^33- (5-14) 

The (incremental) Poisson's ratio 1^12 (and hence z/13 by symmetry) is obtained by setting 
'S'022 = 'S'033 = and using the resulting symmetry L33 = L22 and L22 = —1^12-^11 to obtain 

^^12 = -4oil22/ (-4o2222 + -4o2233)- (5.15) 
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Then, 

5'oil = (^01111 ~ 2z/i2^01122)-^ll (5.16) 

and the (incremental) Young's modulus Ei can be read off as 

El = ^01111 ~ 2z/i2^01122 = -4oilll — 2^01122/ (-^02222 + "4o2233)- (5-17) 

To obtain z/21 = 1^31, 1^23 = z/32 and E2 = -E3, on the other hand, we set 6*011 = 5*033 = 0, 
with Lii = -Z/21L22, -^^33 = -^231,22-, and hence 

^^21-^01111 + ^^23-4-01122 = -^011227 ^^21-4-01122 + ^^23-4-02222 = -4-02233 (5.18) 

5'o22 = (-4o2222 ^ ^21-4oil22 — ^23-4o2233) -^22 = -E'2-^22- (5.19) 



and 
These yield 

and 



_ -4oil22(-4o2222 — -4o2233) _ -4oilll-4o2233 ^ -4oil22 /p- r,n^ 

~ ~A ~A 12 ' ^^23 - J ^2 ^^-^^f 

v^01111vA02222 — -^01122 -^01111-^02222 " -^01122 



-4oilll(-4o2222 + -4o2233) 2^, 



2 

TP _ I A A \ - -Uiiiiv- 'uzzzz 1 - -uzzjjy — ^'01122 (r oi A 

-£^2 — 1-^02222 — -^02233 j 1 1 T2 • l-^./ij 

-^01111-^02222 — -^01122 

The connection 

^2/1^21 = ^i/^^i2 (5.22) 

then follows, as in the classical linear theory. We note, however, that if the increments of 
the Cauchy stress were used in the definitions of the stress-dependent Poisson's ratios and 
Young's moduli instead of the nominal stress (which is entirely legitimate), then this would 
not follow. There are now four independent material parameters: Ui2 = 1^13, 1^21 = 1^31, 
1^23 = ^^32 and El, for example, with E2 = E3 given by fl5.22p . 

For the incremental shear response in the plane of symmetry, we have 

Cr23 = 5023 = -4o2323-^32 + -4o2332-^23 = Q;i(-^23 + -^32); (5.23) 

according to (14.121) 9. and hence ai is the shear modulus in the plane of symmetry. In fact, 
it is straightforward to show that it can be expressed in terms of the other parameters as 

ai = E2/2{l + U23), (5.24) 

similarly to the situation in the classical theory. 

However, when it comes to shear in a plane that contains the axis ei the situation 
differs from the classical one because of the influence of r. We have 

CTi2 = 5*012 = 5*021 + rL2i = [ai + t{(3i + t'Ji) + r]L2i + [ai + r(/3i + r7i)]Li2. (5.25) 

For shear in the xi direction with incremental simple shear Lij = Li2Sii62j we obtain 

cri2 = 5*021 = [ai + r(/3i + r7i)]Li2, (5.26) 
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while for shear transverse to the Xi direction with incremental simple shear Ljj = L2iS2iSij 
we have 

CTi2 = Sou = [ai + r(/3i + r-fi) + t]L2u (5.27) 
and the associated shear moduli are ai + t{/3i + r7i) and ai + r(/3i + rji) + r, respectively. 

Thus, in total, there are six coefficients that are functions of r, but when linearized in 
r there remain six constants, namely ai(0), a2(0), /3i(0), /32(0), a'i(O) ct2(0)5 together 
with T. 

5.3 Planar shear initial stress 

If the initial stress lies in the {xi,X2) plane and is a pure shear stress of amount r then 
r = r(ei (g)e2 + e2 ® ei) and the principal values of r are ±r. The principal axes of r bisect 
the background axes ei and 62, along 612 = (ei ± e2)/-\/2, say. Here, we therefore take as 
our axes of reference the principal axes ei, 62, and §3 = 63; the associated components of 
Aq are 

Aonn = 2ai + ^2 + (1 + 4/3i + 2/32)r + {P^ + 471 + 2'y2y + 2^^t^ + 74r^ (5.28) 
A2222 = 2ai + 02 - (1 + 4/3i + 2/32)r + (/33 + 471 + 2-f2y - 2731-^ + 74r^ (5.29) 
A3333 = 2ai + a2, ^01122 = «2 - f^sT^ + 2-f2T^ + 74r^, (5.30) 
A1133 = "2 + (^2^ + 72^^^ A2233 = "2 - /32T^ + 72^^^ (5.31) 
A1221 = "1 + 27ir^ A01212 = ai + r + 2'jiT^ , A2121 = "i - ^ + 27ir^ (5.32) 
A3131 = A3113 = "1 + + i"^7i, A01313 = ai + r + tI3i + T^ji, (5.33) 
A3232 = A3223 = "1 - + 1"^7l, -4o2323 = tti - T" - + i"^7i- (5.34) 

Components referred to principal axes are indicated by a superposed hat. We also have 
trr = tr(T^) = 0, tr (r^) = 2r^, and in general all the coefficients ai, . . . , 74 depend on r. 

To illustrate the results in this case we consider the restriction to incremental plane 
strain with L^i = Li^ = 0, i = 1,2,3. The components of the incremental nominal stress 
are then given by 

'S'oil = ^01111-^11 + -4.01122-^22; 'S'022 = -4.01122-^11 + -4.02222-^22 ( 5.35) 

for biaxial deformation parallel to the principal axes, and 

5*012 = («i + 27ir^)(Li2 + L21) + rL2i, 5*021 = (ai + 27it^)(Li2 + L21) - rLi2, (5.36) 

for shearing deformations. 

The plane strain Poisson's ratios and Young's modulus Ei are then deduced as 

^^12 = -4oil22/-4o22225 ^^21 = -4oil22/-4oilll , Ei = -4oilll — -4q^]^22/-4o2222 (5.37) 

and, as in fl5.22p . E2 = -E'iz/2i/z/i2. The shear moduh are ai — r + 27ir^ and ai + r + 271 
for shear in the ei and 62 directions, respectively. 

When r = we recover the classical results for plane strain isotropy: Poisson's ratio is 
^^12 = ^^21 = Ao/(Ao + 2/io), and Young's modulus is Ei = E2 = 4/io(Ao + Aio)/(Ao + 2;Uo). 
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6 The effect of initial stress on infinitesimal wave prop- 
agation 



The theory of small amplitude (incremental) deformations or motions superimposed on a 
static finite deformation is well established, but has received relatively little attention in the 
case of an initially stressed material with or without an accompanying finite deformation 
except for works based on Biot's theory in the context of linear elasticity. Here we consider 
incremental motions in an infinite homogeneous medium subject to a homogeneous initial 
stress with Aopiqj having the form given in (14.21) and various specializations of that form. 
From (I3.12P we recall that the equation of incremental motion is 

Consider a homogeneous plane wave of the form 

u = m/(n ■ X — ft), (6.2) 

where m is a fixed unit vector (the polarization vector), / is a function of the argument 
n ■ X — f t with appropriate regularity, n is a unit vector in the direction of propagation, 
and V is the wave speed. Substitution into the equation of motion (16. ip (after dropping 
/", which is assumed to be non-zero) leads to 

AopigjUpngTrij = pv'^rrii. (6.3) 

The associated acoustic tensor Q(n) has components defined by 

Qij{n) = Aopiq-jUpUq (6.4) 

and enables the propagation condition (16. Sp to be written compactly as 

Q(n)m = pt>^m. (6.5) 

For any given direction of propagation n we have a three-dimensional symmetric alge- 
braic eigenvalue problem for determining pv"^ and m. Because of the symmetry there are 
three mutually orthogonal eigenvectors m corresponding to the directions of displacement 
and the (three) values of pv'^ are obtained from the characteristic equation 

det[Q(n) -pt;^I] = 0. (6.6) 

If m is known then pv"^ is given by 

pv^ = [Q(n)m] ■ m, (6.7) 

and corresponds to a real wave speed provided pv"^ > 0, which is guaranteed if the strong 
ellipticity condition holds, i.e. if 

[Q(n)m] ■ m = Aopigjnpnqmimj > for all non-zero vectors m, n. (6.8) 
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Then, a triad of waves with mutually orthogonal polarizations can propagate for any di- 
rection of propagation n. Henceforth we assume that t he strong elliptici t y con ditio n holds. 
For de tailed discussion of strong ellipticity we refer to iTruesdell &: Noll! f jl965l ) and lOgden 
f ll984j ). for example. 

In the following we examine the effect of initial stress on the propagation of plane waves 
for the three examples of initial stress considered in Sections S] and El and for this purpose 
we give explicit expressions for Q(n) in each case. 



6.1 Isotropic initial stress 



For the form of ^Aq given by (14. Sp . with the connections (15. 4p and (15. 5p and (16. 4p . we obtain 
simply 



Q(n) = ^(r)I + [A(r) + yu(r)]n ® n. 



(6.9) 



As in the classical theory of isotropic elasticity with no initial stress there exists a longi- 
tudinal wave with speed vl, say, and two transverse waves with equal speeds vt, say, for 
any direction of propagation. With dependence on r these are given by 



pvl = A(r) + 2/i(r), pv^ = /i(r). 



For 



sufficiently small initial stress we may linearize these expressions to give 

A(r) + 2/i(r) ~ Aq + + [2a'i(0) + ^'2(0) + 4/3i(0) + 2/32(0) + l]r, 
/i(r) ~/io + [«'i(0) + 2/3i(0) + l]r. 



(6.10) 

(6.11) 
(6.12) 



where ai(0) = /iq, ay jO) = An. 



Note that if, as in lTangI fll967| l. the isotropic moduli were defined based on the increment 
in the second Piola-Kirchhoff stress according to (15. Sp then /i(r) and A(r) + 2/i(r) would 
have to be replaced by /x( r) + t and A(r) + 2/i(r) + r, respectively, leading to the erroneous 
conclusion of iTangI ( 19671 ) that when the elastic moduli are independent of the initial stress 
the wave speeds are given by pv^ = fiQ — P, pvl = Aq -l- 2/io — P for the case of a hydrostatic 
pressure (r = —P). Results such as these are not supported by the data shown in Tables 
□ and El 



For a general isotropic elastic material under hydrostatic pressure iTruesdelll (119611) ob- 
tained expressions for the speeds of longitudinal and transverse waves; see also ITruesdell fc Noll 
(119651 ) ■ section 75. 



6.2 Uniaxial initial stress 

For uniaxial initial stress the acoustic tensor is given by 

Q(n) = 5a (g) a + Cn (g) n + ^(n (g) a + a (g) n) 



(6.13) 
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where 



A = ai + [l + (3i + 7iT]r(n ■ a.)^, 

B = PiT + 7ir2 ^ ^ 273r + 74rV'(n " a)', 

C = ai + a2, 

D = [/3i + /32 + (71 + 72)r]r(n ■ a). 



(6.14) 



If n = a (propagation along the direction of unaxial stress), then there exists a longi- 
tudinal wave with speed vu given by 

pvl, = A+B+C+2D = 2al+a2 + (l+4/3l+2/32)r+(/33+47l+272)r2+273r=^+74r^ (6.15) 

and two transverse waves with equal speeds f 12 given by 

pvl^ = A = ai + {1 + l3i)T + -iiT^. (6.16) 



These formulas are consistent with the formulas (5.20) and (5.23) in IShams et al\ ( 1201 ll ) 
relating to propagation along a principal axis of r except that in (5.20) there is a misprint 
(the coefficient of 71 should be 4 instead of 3 — this arises from the fact that in the 
expression (5.18) for a in the latter paper 71 should be 271). The linearized forms of f l6.15p 
and f l6.16p are 



pvl^ = Ao + 2/io + [2a; (0) + ^'2(0) + 4/3i(0) + 2/32(0) + l]r, 



(6.17) 
(6.18) 



respectively. 

On the other hand, if n ■ a = (propagation transverse to the direction of uniaxial 
stress) then 

A = ai, B = (3iT + -fiT^, C = ai + a2, D = 0, (6.19) 

and a longitudinal wave exists with speed f 22 given by pf I2 = A + C. There are also two 
transverse waves, with polarizations along and perpendicular to a and speeds f 21 and f 23 
given by 



A + B = ai + PiT + 7ir , PV23 = A = ai, 



respectively. We notice immediately that 



P^?2 - pvli = r 



(6.20) 



(6.21) 



exactly, as expected from the propertv f|3.8D . This well-known re lationship (see, for ex- 
ample, iBiotl . Il965l : Man &: Lul . 119871 : iHogerl . Il993al : iNorrisl . Il998h forms the basis of an 
experimental acoustic identification of a solid with anisotropy due to an initial stress, 
as opposed to a general anisotropic linearly elastic solid without initial stress, for which 
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PV21 = Cqq (Voigt notation). We note that, even earlier. iBiotI (jl940af ) obtained sepa- 
rate expressions for, in the present notation, vu and f 21 from which the above relationship 
may also be deduced. Finally, we note that the linearized forms of PV22 and fl6.20p are 



PVI2 = Ao + 2po + [2a;(0) + a^(0)]r. 



PV21 



/io + [a'i(O) + 



po + a[{0)T. 
(6.22) 

If the propagation takes place neither along the direction of uniaxial stress nor perpen- 
dicular to it, then several possibilities arise, which are detailed in Appendix A. 



6.3 Planar shear initial stress 

Here we consider wave propagation in the plane of shear — the (xi, X2) plane. The relevant 
components of the acoustic tensor are then given by 

Qu = ai + [ai + as + 2(271 + + 74^>i + 2(1 + 2/3i + (32 + 73i-')™in2 + (/^s + ^-f^ynl, 

Q22 = ai + [ai + a2 + 2(271 + 72)r2 ^ ^^^4]^2 ^ 2(1 + 2/3i + /32 + 73i-')™in2 + + '2li)r^n\, 

Q12 = (2/3i + (32)t + 73r=^ + [ai + 02 + (271 + 272 + /S^y + 74r>in2, 

referred to background axes (not the principal axes considered in Section 133]) . As already 
noted, we assume that the strong ellipticity condition holds so the wave speeds are real. 
As in the previous section we work in the {xi,X2) plane with 77,3 = ma = and set 
rii = cos6,n2 = sin 6' and mi = cos0, m2 = sin(f). Then, by eliminating the wave speed 
from the propagation condition, we may apply equation flA.SP from Appendix A, recast as 

'''' [ai + a2 + (271 + 272 - /33 + 74r2)r2] cos 2^ ' I • J 

which gives for any given 6. 

The (in-plane) wave speeds are given by 



pv 



.2 1 



Qn + Q22±\ {Qn-Q22y + ^Q 
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(6.24) 



For definiteness it is interesting to consider the situation in which r is small and fl6.23p is 
linearized in r, which leads to 

(ai + 02) sin(20 - 29) = 2{2(3i + (32)t cos 20, (6.25) 

from which we deduce that a longitudinal wave can propagate for r 7^ if either 9 = (j) = 
7r/4 or 2/3i + /32 = 0. The first of these possibilities corresponds to propagation along a 
principal axis and the second to a special set of values of the material parameters that 
allows propagation of a longitudinal and transverse wave in any in-plane direction. 
In the linear specialization we obtain 

Qii + <522 = 3ai + a2 + 2(1 + 2^2 + /32)r sin 29, (6.26) 
(Qii - Q22)' + 4Q?2 = ("1 + "2)' + 4(«i + «2)(2/3i + /32)rsin2^, (6.27) 
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and the wave speeds fl6.24p are then given by 



pv^ = Ao + 2/io + {2a;(0) + a^(0) + [1 + 4/3i(0) + 2/32(0)] sin 2^}r, (6.28) 
pt,2 = po + [a;(0) + sin2^]r, (6.29) 

in which the coefficients have now be hnearized in r. These are respectively longitudinal 
and transverse when = 7r/4, as indicated above, or in the special case 2/3i + = the 
tra nsverse wave speed stands but the longitudinal wave speed specializes accordingly. 



Tang fll967l ) considered in-plane wave propagation for an initial shear stress in which 



the only non-zero components of the second Piola-Kirchhoff stress were T12 = T21. As with 
the case of hydrostatic pressure discussed at the end of Section 16. 1[ Tang used an incorrect 
incremental form of the constitutive law. When linearized in the initial stress the results 
from equation (3.10) in his paper that parallel fl6.28p and fl6.29p can be shown to reduce, 
in the present notation, to 

pv^ = Ao + 2/io + (Ao + 3/io)r sin 29, pv^ = /iQ, (6.30) 

where we have set f = t/pq and r = T12 since there is no distinction between stress mea- 
sures themselves in the reference configuration, which is quite different from the situation 
for their increments. Note, in particular, that fl6.30p 9 depends on neither the initial shear 
stress nor the direction of propagation, which is quite unrealistic and cannot be recovered 
from fl6.29p for r 7^ 0. Fortuitously, fl6.30P i can be recovered from (16.280 by making the 
special choices 2a'^(0) + ^2(0) = and 2/3i(0) + /32(0) = 1 + \q/2pq of the coefficients. 



6.4 Some connections with the classical theory of Biot 



It is interesting now to consider how the c l assical theory of in itial stress in the general linear 
theory of elasticity due to Biot (see iBiotl . Il939l . Il940al . Il965l ) relates to the present theory. 
As noted in Section [H the elasticity tensor of Biot, with components Bpiqj depends in an 
unspecified way on th e initially stressed configuration. First, we record that, as shown in 



Qgden fc Singhl (120111 ). the general connection between Aopiqj and Bpiqj may be written in 



the form 



Opiqj 



piqj 



2 "PJ <?« 2 PI *J 



lA .7- . + lA..,- 



8 T • 



(6.31) 



which can be shown to be equivalent to the expression (4.25) given in Chapter 2 of Biot's 
book (Biot, [l965i ). For the general expression (14. 2 p to reduce to the Biot form for isotropic 
response, as quantified in (II. 3p . the material parameters in (14. 2 p must be specialized to 



tti = /io, 02 = Ao, /3i = -1/2, /32 



0, 



(6.32) 



where pq and Aq are the classical Lame moduli and terms of order higher than 1 in r are 
neglected. 

For the hydrostatic stress considered in Section 6.1 these specializations yield, from 
(I6.12P and (16.111) . pv'^ = po and pvl = Aq + 2pQ — r. The first of these agrees with 
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the original result of Biot, who mentioned that any dependence on the initial stress must 
be through the elastic constants themselves. The results in the present paper make the 
dependence explicit. It does not appear that the result pvl = Aq + — t was given by 
Biot. 



We remark that in IShams et al\ ( 120111 ) we adopted a slightly different form of the 



isotropic constitutive law from that given in (II. 3p . namely 

I3piqj = fJ'oi^ij^pq + ^qi^pj) + ^O^pi^qj + ^piTqj , (6.33) 

for which, in the list (16.321) . = is replaced by = 1- 

Turning next to the case of uniaxial stress we find from Section 6.2 first that for 
propagation in the direction of initial stress the formulas (I6.17P and (I6.18P reduce to 
pvfi = Ao + 2/io — r and pvf2 = /^o + ^t, respectively, while for propagation normal to the 
initial stress the formulas ( I6.22p reduce to pffg = Aq + 2/io, P'^^ii = /^o — |t and pffg = po- 
For the case of compressive initial stress with r = — P (P > 0) the speeds of the relevant 
two transverse waves agree with those obtained by Biot, specifically pv^ = po± |P. 

It is also interesting to apply the Biot specialization to the shear stress example in 
Section 6.3. Equations (E28D and ( K29\\ yield 

pvl = Ao + 2/io — T sin 26*, pv^ = po + t sin 29, (6.34) 

and a longitudinal and transverse wave can propagate in any in-plane direction, where 6 
is the angle the propagation direction makes with the principal direction of stress corre- 
sponding to principal stress +r. Note that for propagation along either principal direction 
there is no influence of r, which would seem to be unrealistic. 



7 Deformed and pre-stressed isotropic elastic solid 

In this section we consider the initial stress to be associated with a finite deformation from 
an unstressed configuration (denoted Bq) of an isotropic elastic solid, and we shall con- 
sider two states of the accompanying stress analogous to those considered in the foregoing 
sections, specifically pure dilatation, corresponding to hydrostatic stress, and an axially 
symmetric deformation corresponding to uniaxial stress. There is no direct analogue of the 
planar shear stress situation since when accompanied by deformation such as simple shear 
there will in general also be normal stresses, which are not considered in Section 14.31 

For an isotropic elastic solid the Cauchy stress tensor cr is given by the appropriate 
specialization of (12.141) . namely 

Ja- = 2WiB + 21^2 (/iB - B^) + 2/31^3!, (7.1) 

where we recall that J = det F and B is the left Cauchy-Green tensor, where F and J are 
now measured relative to Bq instead of Br- We may consider the strain energy to depend 
on the three principal stretches Ai,A2,A3 instead of the principal invariants Ii, 12,13, and 
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for this purpose we write W = W{\i, A2, A3). The principal Cauchy stresses cxi, CT2, (J3 are 
then given simply by 



JcTl = XlWi, J(T2 = X2W2, J(T3 



(7.2) 



which are equivalent to the principal components of (17.11) . where Wi = dW/dXi, i = 1,2, 3, 
and J = A1A2A3. This is easily seen by noting that in terms of the principal stretches the 
invariants /i, /2, -^3 are given by 



h — \l + \l + A3, I2 — A2A3 + A3A^ + A^A2, I3 — — A^A2A3. 
In terms of invariants the components of the elasticity tensor are given by 



(7.3) 



JA 



Opiqj 



2{Wi + IiW2)Bpg5ij + 2W2[2BpiBgj — BigBjp — BprBrgSij — BpgB^j^ 



+ 2l3W3{2SipSjg - Sig6jp) + 4WuBipB,g + 4W22{hBip - Bi,Brp){hB,g - 
+ 4Wi2{2IiBipBjg — BipBjrBrg — BjgBirBrp) + AIsWis{Bip6jg + Bjgdip) 



BjgB gg) 



+ 4.hW2z[h{Bip5jg + Bjgdip) - 5ipB,jr.Brg - 6jgBirBrp] + il^WssSir^Si 



(7.4) 



where Bij are the cor aponents of B. Thi s specializes the form of JAopigj for an initially- 
stressed solid given in IShams et al\ (1201 ll ) to the present situation, but equivalent forms of 
fl7.4D can be found in the earlier literature on isotropic finite elasticity (see, for example. 



Haves fc Rivlinl . Il96lh . 



When referred to the principal axes of B the components (17. 4p can be expressed more 
compactly as 



XiWi 



■Oijij 



2 

A2-A| ^' 



Oijij 



(7.5) 



where Wij = d'^W /dXidXj, i,j G {1,2,3}. Note that the only non-zero components of 
Afi are Annji , i,j G {1,2,3}, together with, for i 7^ j, Aoijij and Aoijji (see, for example. 



Qgdenl . Il984[ ). as is the case for the components of the elasticity tensor of an initially 



stressed but undeformed material when referred to the principal axes of the initial stress 
T, as can be seen by reference to (14. 4p . 

If two of the principal stretches coincide, say Aj = Aj, then aj = ai and a limiting 
process can be applied to express the second and third elements in (17. 5p as 



Aoijij y 2 iAoiiii At 



■Oiijj ~l~ ^i)' Aoijji y ^^i^AQiiii Aoiijj Oi 



(7.6) 



respectively. 



7.1 Pure dilatation 

Here we consider a pure dilatation with corresponding isotropic stress. Let Aj = J^^^, i = 
1,2,3, and define W{J) = W{J^^^, J^^^, J^^^). We denote the corresponding (equal) prin- 
cipal stresses by a . Then the simple connection a = W'{J) follows. The only independent 
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components of ^ are then Aoujj = J ^^^Wij evaluated for the pure dilatation, with Aouu 
independent of i and Aoujj independent of i and j 7^ i. Now, from flS.ip . we obtain 

& = AqL + (tL - (T(trL)I. (7.7) 

We consider an incremental simple shear deformation, and, without loss of generality, we 
may take this to be in the (xi, X2) plane. Then 

1 

0"12 = -^01212 (-^^21 + L12) = -{Aouu — Aonjj + Cr)(-Z^21 + L12) . (7-8) 

This enables us to define the incremental shear modulus as a function of J, which we write 
as /i(J). It is given by 

KJ) = 2(^1111 - A1122 + (t). (7.9) 

Note that 6*021 = ^02121-^12 for simple shear in the Xi direction, 5912 = ^01212-^21 for 
simple shear in the X2 direction and ^02121 = "4oi2i2- We note in passing that T012 = 
(-^01212 — cr)(-^2i + L12) and it would be incorrect to define the shear modulus as ^01212 — cr 
based on use of To- 

Next consider an incremental pure dilation e, so that Ln = L22 = -Z^ss = e/3 and 

tro- = (Ann + 2An22 - 2a)e. (7.10) 

The incremental bulk modulus, which we denote by k{J), may be defined as a function of 
J as 

k{J) = -(Ann + 2An22 - 2cx), (7.11) 

where the components of Aq are again evaluated for Aj = J^^^, i = 1,2, 3. This may also 
be expressed in the form k (J) = JW"( J), which agrees with a definition of the incremental 
bulk modulus adopted by IScottI (120071 ). 

If we set a = T, where r is the initial hydrostatic stress considered previously then, 
because of the connection r = W'{J), we can in principle switch between the two different 
formulations, although, for a given r, this would involve inversion of the relation r = W'{J) 
to obtain J. When the switch is made we can identify fi{J) and k{J) with /i(r) and 
k{t), respectively. We should note here that the densities in the stress- free reference 
configuration Bq, with density po, and the deformed (or initially-stressed) configuration 
B = Br are related by po = pJ and the factor J needs to be used to switch between pov"^ 
and pv"^ in considering formulas for various wave speeds, where p = p^.. 



24 



7.2 Uniaxial stretch with lateral contraction 



Now consider a uniaxial stress ai = a with a2 = cs = and stretches Ai and, by symmetry, 
A2 = A3. Then, W2 = W3 = and the components of are given by 



^"^01111 
JAoiui '- 



AiWll, ^^02222 — ^\W22 



JA 



03333, 



AiA2iyi2, JA 



Olili 



aFaI' 



JA 



■Oilil 



JAq223?, — A 2^^23 5 
AiA^W^i 



A? - A2 



2,3, 



Aquh — Aquu - 
v4.02323 = "4.03232 



Aquh — Aquu — 0", 



2,3, 



"^02332 — "^03223 — -^{Aq2222 — "^02233) 



(7.12) 
(7.13) 
(7.14) 
(7.15) 



all evaluated for A3 = A2 and with o = A2^iyi. Poisson's ratios and Young's moduli can 
be defined in exactly the same way as in Section 5.2, except that here the components of 
Aq are different from those in Section 5.2. There is no need to repeat them all here, but 
we note, for example, that 



^12 — '4oil22/(-4.02222 + -^02233) — AiVri2/A2(W22 + ^^23^ 



Ei = A 



01111 



2Z/12A1122 = J-'Xl[Wu - 2W^J{W22 + W23)], 



(7.16) 
(7.17) 



which are both functions of Ai when A2 = A3 is determined from W2 = for a given form 
of W. The ex pression for 2 is consistent with the definition of the incremental Poisson's 
ratio given by IScottI (120071 ) . On the oth er hand, the definition of Ei above differs from the 
corresponding definition in IScottI ( 120071 ) since the latter is defined, in the present notation, 
as Aidcr/dAi with A2 = A3 given by W2 = 0. The definition of Ei above corresponds 
to AidH^i/dAi updated to the deformed configuration by the push-forward factor J~^Ai, 
which is the appropriate specialization of the general push forward operation J~^F that 
takes the nominal stress S to the Cauchy stress cr = J^^FS. This definition of Ei based 
on nominal stress is consistent with that used in Section 5.2. 



7.3 Application to second-order elasticity 

For definiteness we now specialize the form of strain-energy function W to the third order in 
the strain. The precise form of this approximation depends on the choice of strain measure, 
but here, in order to make contact with several contributions to the literature, we shall use 
the Green strain tensor E = |(C— I), where C is again the right Cauchy-Green deformation 
tensor. The aim here is to obtain the first order correction to the classical longitudinal 
and transverse wave speeds. For a discussion of the advantages of using logarithmic strain 
instead of Green strai n, particularly when appro aching the incompressible limit, we refer 
to the recent paper bv lDestrade fc Qgdenl (l2010f ). 
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7.3.1 Historical expansions of the strain energy 

The use of inv ariants for expre ssing the third-order expan sion of W appe ars to have been 
introduced by iBrillouinI ( 119251 ): see also his monograph f lBrillouinl . Il946l ). The Brillouin 
expansion may be written in the form 



W = Wo- Ipoh + + 7,^^0/2 + AlJ, + Bil + Ch, 

Z o 4 



(7.18) 



where, in different notation from that used by Brillouin, Ji = 2trE, I2 = 4tr(E^), I3 = 
8tr(E^), respectively of orders 1, 2, 3 in E, the constant Wq is the energy in the reference 
configuration and po corresponds to an initial pressure in the reference configuration. Note, 
in particular, that Brillouin used 2E rather than E itself as the strain measure and that 
the notations A, B, C are different from those defined in (16.141) . 

Let us drop the terms Wq and po, which are redundant for our purposes, and recast 
the remaining terms usin g the principal invar i ants o f Green strain, which, for consistency 
with the notation used in lDestrade fc Ogderu (j201Cll ). we denote by 21,^2, "^3- Thus, 



trE, 



is = det E. 



(7.19) 



Then, we have 



W = -(Ao + 2/io)i? - 2fioi2 + 8iA + B + C)il - 8{2A + 3C)iii2 + 24Ci3. (7.20) 

This is entirely equivalent to the strain-ener gy function genera l ly ref erred to as the Mur- 
naghan form of strain energy that appears in iMurnagharu (jl95ll . 119671 ) and is based on the 
use of Green strain. It is written 



W = ^{Xq + 2fio)il - 2/io«2 + + 2m)i\ - 2miii2 + ni^, 



(7.21) 



where l,m,n are the Murnaghan constants. The Brillouin constants A,B,C are given in 
terms of /, m, n by 



A = —(2m-n), B = —(21 - 2m + n), C = —n. 
16^ ^' 48^ ^' 24 



(7.22) 



We remark that in his original paper iMurnaghan fll937h worked in terms of the Almansi 
strain tensor (I — B~^)/2 and its principal invariants, which we denote here by Ji, J2, J3. 
The original energy function of Murnaghan has the form 



W = -(Ao + 2^o)/i - 2/^0/2 + 11! + rfihh + nh. 



(7.23) 



In general (17.231) is different from fl7.2ip . but the two are equivalent to the third order in 
the strains. It is then easy to show that the constants /, m, n and Z, m, n are related by 



r=-(/ + 2Ao) 



m 



-2m - 4Ao - 12/io 



n 



n + 12/io 



(7.24) 
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This difference has significance when considering approximations to wave speeds at this 
order, as we shall see shortly. We note in passing that because of the requirement of 
objectivity the Almansi strain tensor (or any other Eulerian strain tensor) can be used as 
the argument of the strain energy function if and only if the material is isotro pic. Clearly, 
the wo r k of Brillouin on this top ic has been to sq r ne ext ent overlooked , although iMurnaghan 



fll937n . iHughes fc Kellvl fll953f ). and iTruesdelll fll96lh did refer to iBrillouinI (119251). For 



furthe r historical discussion of s econd-order elasti c ity, in cluding the contribution of iRivlin 
(Il953l ). we refer to section 66 of iTruesdell fc Nolll (Il965[ ). 

An equivalent for r n of t he third-order expanded energy function was also introduced 
by lLandau fc Rumerl (jl937l ). who were apparently unaware of the work of Brillouin. This 
may be written 



W = ^Ao(trE)2 + /iotr(E2) + \Aii{Y?) + S(trE)tr(E2) + ^^(trE)^, 



(7.25) 



where o yerbars have been u s ed to distinguish the third-order constants from those of Bril- 
louin; in iLandau fc Lifshitzl (119861 ) the notations A, B, C were used, differing from those in 
Landau fc Rume'ir ( 1937l )~~The connections between the Murnaghan c onstants l,Tn,n and 
the Landau constants A, B, C were noted in lDestrade fc Ogderu (120101 ) as 



A 



n. 



B 



m n, 

2 ' 



C 



I — m -\ — n. 
2 



(7.26) 



BiotI ( ll940bl ) also developed a third-orde r expansion , which is equ ivalent to the above 



and details may also be found in his book ( iBiotl . Il965[ ). iBiotI (Il965[ ) worked in terms of 
the principal strain components Aj — 1 and the corresponding principal Biot stresses. His 
third-order constants, denoted D, F, G, can be shown to be related to /, m, n via 



= / + 2m + -(Ao + 2/io), F = / + -Ao, G = 21 - 2m + n. 



(7.27) 



He did not give the form of strain energy explicitly. 

Finally, we mention the third-order expansion adopted by lToupin fc BernsteinI (jl96l[ ). 
who used the invariants trE, tr(E^), tr(E^) and third-order constants Ui, i>2, z/3. In terms 
of the principal invariants of Green strain their energy function has the form 



W 



1 



(Ao + 2fio)il - 2/ioi2 + ^('^i + 6z/2 + 8z/3)if 

6 



2(z/2 + 2z/3)zii2 + 4z/3i3. 



(7.28) 



It i s straightforw ard to show that the constants z/i, z/2, z/3 are related to /, m, n and A, B, C 
by (iNorrisl . [l998h 



ui = 21 



2m + n = 2C, = m-\n = B, u^ = \n = \a. (7.29) 



2 '-^44 

Connections betwe en some of the ab o ve set s of constants and others used by RivlinI ( 1953 ) 
were also noted by ITruesdell fc Nolll (119651 ). 
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Here we shall work in terms of the Murnaghan constants l,m,n but cast ( I7.2ip in terms 
of the principal invariants fl2.9p of C as 

W = ^{h-3r + ^{ll-2h-2h + 3) 

o 4 

/ 777 77 

+ ^(/i - 3)=^ + Y^(A - - 3h) + -ih -h + h- 1). (7.30) 

For this energy function we have W22 = W^is = W^23 = W33 = and the expression f l7.4p 
reduces to 

JAopiqj = 2(Wi + IiW2)Bpg5ij + 2W2[2BpiBqj — BiqBjp — BprBrq5ij — BpgBij] 

+ 2hWsi26,p6,q - 6,q6,p) + AWnB.pB^q 

+ 4:Wi2{2IiBipBjq — BipBjrBrq — BjgBirBrp), (7-31) 

and the remaining coefficients Wi, W2, W3, Wu and W12 are simply obtained from fl7.30p . 

In working with second-order elasticity the corrections to the classical moduli are ob- 
tained by expanding the coefficients in the above to the first order in E. We have C = I+2E 
and Ii = 3 + 2E, exactly, which we use together with the linear approximations I2 — 3+4:E, 
Jg ~ 1 + 2E, B ~ C = I + 2E, where E = trE. We also note that p ~ pr{l - E). To the 
first order in E we obtain 

W^i = /io + |n+ i(Ao + 2/io + 2m)E, W2 = -^Po - - ^mE, W3 = ^n, 
Wn = \iXo + 2/io + 4m) + i(/ + 2m)E, W,2 = -\m, (7.32) 
and hence 

J-Aopiqj ~ PoiSijSpq + 6iq6jp) + Xo^ipSjq + ^ [(2Ao + 2m - n)6ij6pq + 2(2/ - 2m + n)6ip6jq 
+ (2m - n)Siq6jp]E + ^(4/io + n){6pqEij + SiqEjp + SjpEiq) 

+ ^(8/^0 + n)6ijEpq + (2Ao + 2m - n){6ipEjq + 6jqEip). (7.33) 

We are particularly interested in the case of a pure dilatation, for which, with Eij = 
E6ij/3, (17.331) reduces to 

JAopiqj = Po{Sij6pq + Siq6jp) + XoSipSjq + ^(12/io + 6Ao + 6m - n)E6ijSpq 

+ ^(8/io + 6m - n)E6iq6pj + ^(4Ao + 6/ - 2m + n)E6pi6qj. (7.34) 
From the definitions fl7.9p and (17. lip we now obtain the approximations 

/i(J) ~/io + 7(6Ao + 6/io + 6m-n)^ (7.35) 
6 
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and 

k{J) Ko + -{9l + n)E, (7.36) 
9 

for small dilatation, where to the first order J = 1 + E, and we note that fi{l) = /xq, 
k{l) = Kq. 

To the second order we may expand the isotropic stress as 

a = W'{J) ~ eW"{l) + ^e^W'"{l), (7.37) 

where W"{1) = kq, e = J — 1 — E + and, for the Murnaghan strain energy, W"'{1) = 
-Ko + 21 + 2n/9. 

As indicated earlier, if the initial stress r discussed in Sections H] and O is associated 
with a pure dilatation, so that r = rl, then in principle the results here can be converted 
to those based on r. In particular, if we set a = r in the above and work to second order 
then we may invert the r < — t- e relation in the form 

e ^ t/ko - ^W"\iy/Kl (7.38) 

but in considering the linear approximations of the shear and bulk moduli only the first 
order term need be retained. Then, with r = kqE, we can identify /u(r) and k,{t) with fi{J) 
and k(J), respectively. Thus 

/x(r) ~ /io + [ai(0) + 2/3i(0) + l]/to£ ~ fi{J) ^ /io + ^(6Ao + 6/io + 6m - n)e, 

6 

K(r) ~ Ko + ^[2a'i(0) + 4/3i(0) + 3a^(0) + 6/^2(0) - l]Koe ^ k{J) ~ + ^(9/ + n)e. 

(7.39) 

Hence we can relate the constants a;'i(0), 0:2(0), /3i(0), /32(0) to the Murnaghan constants. 
After a little rearrangement, this yields 

a[ (0) + 2/3i (0) = (2/io + 6m - n) /6ko, (7.40) 

a'M + 2/32(0) = (Ao + 6/ - 2m + n)/3fi:o. (7.41) 

In fact, the separate values of a'i(O), 0:2(0), /3i(0) and /32(0) can be obtained by con- 
sidering other initial deformations than pure dilatation, and we show this in Section 7.3.3 
by comparing wave speeds based on the theory of uniaxial initial stress from Section 6.2 
with those based on a finite deformation from an isotropic reference configuration under 
uniaxial stress. 

7.3.2 Implications for the wave speeds 

Turning now to the propagation of plane waves, we note that for any (m, n) pair satisfying 
fl6.5p the wave speed v is given by 

PqV^ = J[Q(n)m] ■ m = J AopigjUpngminii, (7.42) 
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where the factor J is now included to reflect the change in reference configuration from Bq 
to Br = B. For the case of pure dilatation we then obtain, on use of fl7.34l) . 



Po + {Po + Ao)(m • n) 



g(12/io 



6Ao + 6m — n)E 



+ -[8Ao + 8/xo 



For a longitudinal wave with m = n this reduces to 

1 



12/ + 2m + n\E{m ■ n) 



Povl 



Ao + 2/io 

and for a transverse wave with m ■ n = 

Po + (Ao + 2pQ + m 



-(7Ao + 10/io + 6/ + 4m)E, 



PqVt 



'-n)E. 



(7.43) 



(7.44) 



(7.45) 



These results agree with those obtained by iHughes fc Kellvl fll953l ) for the case of hydro- 
static pressure (r = — P). They used th e Murnaghan energy function based on Green 



strain. As shown by IShams et al\ ( 20111) there is s i milar a greement in the case of the 
uniaxial compression considered by IHughes &: Kellyl (119531 ). iToupin &: BernsteinI ( Il96ll ) 



obtained equivalent results based on their third-order expansion, but expressed in terms 
of the acoustoelastic coefficients, which we write as 



Je 

d 

AE 



(Pot^i)l 
(Po4) 



E=0 



;(7Ao + lO/iQ + 3z/i + lOz/2 + 8Z/3), 
;(3Ao + 6/io + 3z/2+4z/3). 



(7.46) 
(7.47) 



1 

3' 

_ 1 

Note that the E in Toupin &: Bernstein ( 1961 ) is 3x that used here. Toupin &: BernsteinI 
(|l96ll ) mentioned that their results were equivalent to those obtained by iBrillouinI (|l925l ). 
The Brillouin results therefore pre-date much of the subsequent work. 
In calculating the second order correction to longitudinal and transverse wave speeds. 



Birch! (119381 ) used the original (Almansi strain based) form of the Murnaghan strain energy 
but effectively set I = rh = n = 0, obtaining the results 



Povl 
PoVt 



Ao + 2fio + P(13Ao + 14/io)/(3«:o) = Aq + 2^o - (13Ao + Upo)E/3, (7.48) 
/io + P(Ao + 2i2o)/ko = Po- (Ao + 2^o)^ (7.49) 



for the wave speeds, where P = —hqE is the pressure. The omission of the third-order 
constants can be significant since typically they are of the same order of magnitude as the 
Lame moduli Aq and /^O) as illustrated by the data shown in Table [31 
If the third-order constants are retained then we have, instead, 



Po^l 



PoVt 



Ao + 2/io - (13Ao + 14^0 - 18/ + 2m)E/3, 

Po - (Ao + 2yUo + l-rh + ^n)E. 

2 6 



(7.50) 
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Table 3: Lame constants and Landau third-ord e r elas tic moduli for five solids (expressed 
in units of 10^Nm~^j, as collected by Poruboi i 200d) : his Murnaghan constants m, n, I 
have been converted here to the Toupin-Bernstein constants vi, V2, 



Material 


Ao 






V2 




Polystyrene 


1.71 


-21.2 


-10 


-8.3 


-2.5 


Steel Hecla 37 


-354 


82.1 


-358 


-282 


-89.5 


Aluminium 2S 


-204 


27.6 


-228 


-197 


-57 


Pyrex glass 


264 


27.5 


420 


-118 


105 


Si02 melted 


72 


31.3 


-44 


93 


-11 



If, by contrast, we set the third-order constants /, m, n to zero then the lHughes fc Kelly 



( 119531 ) results reduce to 

P^vl = Ao + 2/io + (7Ao + 10/io)E/3, 

Po4 = /^o + (Ao + 2/io)^. (7.51) 

Note, in particular, the opposite sign but equal magnitude of the second term in the 
shear wave expression compared with f l7.49p . Thus, interpretation of the results requires 
caution. In particular, care must be taken that the results allow for an increase as well as 
for a decrease of the wave speeds with pressure and uniaxial stress, depending on which 
solid is considered (see Tables [1] and [2]) . 



7.3.3 The case of uniaxial stress 

We now consider a deformation from a stress-free configuration of an isotropic material 
associated with a uniaxial stress a in the ei direction. We denote the corresponding 
component E\\ of the Green strain tensor by E. Then, by setting the lateral stress to 
zero and by symmetry, E22 = = — Aoi?/2(Ao -l- /xq), to the first order in E, and 
a = ^iKqUqE / (Ao -l- /io), and hence 

J = 1 + fXoE/{\o + /io) = 1 + ct/Sko, (7.52) 

also to first order. 

The wave speeds f n, fi2, "^22, "^^23 are then calculated by using fl7.33p and appropriate 
specializations of fl7.42p and the connection po = pJ- After some manipulations, which are 
omitted, this yields the formulas 

pt-n = Ao + 2po + 2[2Ao + 7Ao/io + 4po + Po^ + 2(Ao + po)m]cr/3Kofio, (7.53) 

pfi2 = Po + [4po(4Ao + 3po) + 4pom + Aori]o-/12Kopo, (7.54) 

PV22 = Aq + 2po - (2Ag + 5AoPo + 2po - 2po/ + 2Aom)(T/3Kopo, (7.55) 

P^23 = Po- [2po(2Ao + po) - 2pom + (Ao + po)n]cr/6Kopo, (7.56) 
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and we also note the connection pf ^2 ~ — ^- By setting a = r we then compare these 
results with the formulas in (16.171) . (16.181) and (I6.22p . which we now collect together as 



pvl, = Ao + 2/io + [2a; (0) + a^(0) + 4/3i(0) + 2/32(0) + l]r, (7.57) 

PVu = Po + [«i(0) + /3i(0) + l]r, (7.58) 

pvl^ = Ao + 2/io + [2a; (0) + a^(0)]r, (7.59) 

pt4 = /io + «;(0)r, (7.60) 



with PV21 given by (I6.2ip . Note, in particular, that bars have now been placed over a;(0) 
and 02(0). This is because the arguments of ai and 0:2 are different for hydrostatic stress 
and uniaxial stress, respectively the relevant invariants of r are (3r, 3r^, 3r^) and (r, r^, r'^), 
so that a;(0) = 3a; (0) and 02(0) = 3a2(0), while /3i(0) and /32(0) are the same in each 
case. 

Comparison the two sets of formulas yields the results 

/3i(0) = 1 + n/4po, /32(0) = (2Ao + 2m - n)/2po, (7.61) 
"'i(O) = -[2(2Ao + po)po - 2/iom + (Ao + fio)n]/2Kopo, (7.62) 
"2(0) = -[(2Ao + /io)Ao - 2po/ + (2m - n)(Ao + po)]/Kopo, (7.63) 

from which it is easy to check that the results (17.401) and (17.411) are recovered. 

It is interesting that the four constants a[{0), 0:2(0), /3i(0) and /32(0) are expressed in 
terms of the three Murnaghan constants. This is explained by noting that the anisotropic 
constitutive law for an initially stressed material with no accompanying deformation is 
specialized to isotropy by introducing a stress- free reference configuration and an associated 
deformation. 

Finally, we note that the Biot values (16.321) are obtained from the latter formulas by 
specializing the Murnaghan constants to m = —21 = — (Ao + 2/io), ''^ = —Qpo- 



8 Concluding remarks 



In this paper we have examined in detail the effect of initial stress on the propagation of 
small amplitude homogeneous plane waves in an undeformed elastic material on the basis 



of the general theory of a hyperelastic materi al with i nitial stress develope d by lShams et al 



fl201ll ). which had its genesis in the work of iHogerl fll985l . Il986l . Il993al jbl) in particular. A 
key feature of the constitutive law, formulated in terms of invariants of the deformation 
and initial stress, is that the elasticity tensor depends in general in a highly nonlinear way 
on initial stress. Important special cases considered within the general framework included 
initial stresses corresponding to hydrostatic stress, uniaxial stress, and shear stress, for 
which explicit and relatively simple forms of the elasticity tensor were given. For each 
of these states of stress the dependence of various elastic moduli on the initial stress was 
made explicit. For example, simple formulas were obtained for the stress-dependence of 
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the Lame moduli in the case of isotropic initial stress and Poisson's ratios and Young's 
moduli for the cases of uniaxial initial stress and planar initial shear stress. 

The results were applied to infinitesimal wave propagation and it was shown how some 
known results fit within the general framework, and some discrepancies in some of the 
earlier work were highlighted. We then considered the initial stress to be a pre-stress as- 
sociated with the deformation of an isotropic elastic material from a stress-free reference 
configuration and made connections with the analysis from the preceding sections. Specif- 
ically, we considered a pure dilatational deformation and a deformation corresponding to 
simple tension. We then discussed the specialization of second-order elasticity in detail 
a nd collated vario us contributions from the earlier literature that date back to the work 
of iBrillouinI (119251 ). with particular reference to expressions for longitudinal and transverse 
wave speeds, again showing how the results are captured within the general framework 
herein. 
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Appendix A. Non-principal waves in a solid under ini- 
tial uniaxial stress 

Here we complete the analysis of Section l6^ by considering the case of non-principal wave 
propagation, for which the direction of propagation n and the direction of uniaxial stress 
a are neither parallel nor orthogonal. The constants A, B, C, D are as defined in fl6.14p . 
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By solving the equation (16. 6 p with Q(n) given by fl6.13p we find that the wave speeds 
are given by 

pv^ = A, {A-pv^f + [B + C + 2D{n-Si)\{A-pv^) + {BC-D^)[l-{ii-a.f]={). (A.l) 
We therefore consider the cases pv"^ = A and pv^ ^ A separately. 

Case 1: pv'^ = A. The propagation condition (16.51) yields 

[B{m ■ a) + D{m ■ n)]a + [C(m ■ n) + D{m ■ a)]n = 0. (A.2) 

If C(m ■ n) + ^(111 ■ a) 7^ then n = ±a and hence {B + C ± 2D){m ■ a) = 0. The 
case m • a = was covered in Section 16.21 , but there is now an additional possibility, that 
B + C ± 2D = 0. Both these options lead to the same result, which, on substitution from 
(I6.14p . is written 

ai + a2 + (3/3i + 2/32)r + (371 + 272 + /33)r^ + 2-f3r^ + -f^r^ = 0. (A.3) 

There is no restriction on the direction of polarization m. Note that for the specialization 
(I6.32P this yields r = — 2(Ao + po) and A = — Aq and for several of the values of Aq listed 
in Table [3] there is no real wave speed in this case. 

Next, consider the possibility that n 7^ ±a. Then, if C(m ■ n) + D{m ■ a) = it follows 
that also B{m ■ a) + D{m ■ n) = 0. By combining these we deduce that 

{BC - D^){m- a) =0, {BC - D'^){m ■ n) = (A.4) 

provided C ^ 0, D ^ 0. Then, if BC — 7^ we must have m ■ a = and m • n = and 
m is normal to the plane of a and n. Thus, a transverse wave exists for any direction of 
propagation. On the other hand, if BC — = then n is determined from the equation 

BC-D^ = (ai + a2){(3iT + j.r'' + [/Jg + 273r + 74r']r'(n ■ a)^} 

- +P2 + ill + 72)r]V2(n • a)^ = 0. (A.5) 

Since we are considering the case n 7^ ±a and n ■ a 7^ 0, possible directions n generate a 
cone with axis a, provided |n-a| < 1. For each such n, m must satisfy m- {Ba + Dn) = 0. 
Note that in the linear approximation (lA.Sp cannot hold unless /3i(0) = 0, in which case n 
is unrestricted and A = /iq + Q;'i(0)r + r(n ■ a)^. 

Other special cases are as follows: if -B 7^ 0, C 7^ 0, = then m-n = 0, m-a = and 
either n-a = or /3i + /32 + (71 + 72)^ = 0; if 5 = 0, C 7^ 0, -D = then m-n = and either 
n ■ a = or /?! + /32 + (71 + 72)^ = 0. In the latter, if i? = and n ■ a = then (3i + 7ir = 
(assuming, of course, r 7^ 0), while if i? = and n ■ a 7^ then (n ■ a)^ is determined from 
B = 0. Finally, we note that if C = D = then there are four possibilities: (i) n-a = and 
B = and hence /3i + 7ir = — there is no restriction on m; (ii) n - a = and m - a = — 
this is captured by the discussion around (I6.19P : (iii) /3i + /32 + (71 + 72)^ = and 5 = 0, 
the latter determining (n ■ a)^ — there is no restriction on m; (iv) + (32 + (71 + 72)^ = 
and m ■ a = — there is no restriction on n. 
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Case 2: pv'^ ^ A. Now, from the propagation condition f l6.5p with the speciahzation 
fl6.13p we have 

m ■ n)]a + [C(m • n) + D(m ■ a)]n = 0. (A. 6) 

Thus, m, n, a are coplanar unless either the coefficient of a or n vanishes. If neither of the 
coefficients vanish then, without loss of generality, we may confine attention to the (xi, X2) 
plane and set a = ei. Let (ni,?T,2) and {1711,7712) be the in-plane components of n and m, 
respectively and set = = 0. The propagation condition f l6.5p then specializes to the 
two components 

Qumi + Qi2m2 = pv^nii, Q^rrii + Q22"^2 = pv^m2. (A.7) 

For a given propagation direction the wave speed is given by one of the two solutions of 
the quadratic in (lA.ip and is known explicitly. Elimination of ^ from (IA.7P then gives 

{Qu- Q22)mim2 = Qi2{m\- ml), (A. 8) 

which determines the polarization m. Let ni = cos 6, 77-2 = sin 6 and mi = cos 0, m2 = sin (p. 
Then the above can be rewritten as 

C sin 20 + 2D sin ^ , 
= B + Ccos29 + 2Dcos9- ^^-'^ 

From this we can immediately recover some of the previous results. If, for example, 6 = 
(propagation in the direction of initial stress) then (p = or it/ 2, corresponding to 
longitudinal and transverse waves, respectively, with wave speeds given by pvl = A + B + 
C + 2D and pv^ = A (degenerate case). If 6 = 7i/2 (propagation transverse to the initial 
stress) then D = and again = or 7r/2, transverse and longitudinal, respectively, with 
wave speeds given by pv^ = A + B and pvl = A + C. There is, however, an additional 
case not covered previously in which a longitudinal and transverse wave can propagate. By 
setting = 6^ in (lA.Qp and discarding cases already discussed we obtain B cos 9 + D = 0, 
which leads to (on discarding a factor t 7^ 0) 

2/3i + (32 + (271 + 72)r + (/33 + 273r + ^,t^)t cos^ 6 = 0. (A.IO) 

If this has a solution (or solutions) 6 for cos^ < 1 then a longitudinal wave can propagate 
in the direction defined by such an angle (or angles). The corresponding wave speed is 
given by 

pvl = A + C + Dcose = 2ai + a2 + [l + 2/3i + P2 + (271 + 72)r)]r cos^ 6. (A. 11) 

Equally, by setting = 6 + ti/2 the same reduction Bcos9 + = is obtained and a 
transverse wave can accompany the longitudinal wave and has wave speed which, on use 
of (lA.lOp . can be written 

pv^ = A + B + DcosO = «i + /3ir + 7ir2 + rcos2 0-[2/3i + /32 + (27i + 72)r]rsin2 0. (A.12) 
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For the specialization fl6.32p equation (lA.lOp is satisfied and equations (lA.lip and (1A.12P 
reduce to pvl = Aq + 2fiQ + r cos^ 6 and pv^ = /io — |t + r cos^ 9, respectively. Formulas in 
Section 6.4 are recovered by taking ^ = and 9 = ^vr. 

With reference to (IA.6P we now consider the special cases corresponding to vanishing 
of one or other of the coefficients of a and n. These are B{m ■ a) + D{m ■ n) = with 
C(m ■ n) + D{m • a) 7^ and B{m ■ a) + D{m ■ n) 7^ with C(m • n) + D{m • a) = 0. The 
first of these corresponds to the case B cos 9 + D = just considered, while for the second 
m is aligned with the direction of uniaxial initial stress and C cos 9 + D = 0, which yields 
the nontrivial result 

ai + aa + + 1^2 + (71 + 12)t]t = 0. (A.13) 
This puts no restriction on the direction of propagation n and the wave speed is given by 

pv^ = A + B + Dcos9 (A.14) 

as in flA.12p . but the wave is not necessarily transverse in this case. For the special case 
fl6.32p the condition flA.13p yields r = — 2(Ao + Po), as for (IA.3p . and pv'^ speciahzes to 
f'O + (Ao + /io) sin^ 
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